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1. INTRODUCTION 



An essential feature of transonic "axial" compressor stages is the 
contraction of the flow passage approaching and through the blading. 

This is necessary since the stage pressure ratio is significantly greater 
than unity. In order to design blading, or to calculate inlet conditions 
during operation of a compressor, the radial distributions of velocity 
and flow angle are required in a radial plane immediately ahead of the 
blading. Axisymmetry is assumed. 

The purpose of this report is to document calculations performed 
for the "TRANSX" compressor in the Turbopropulsion Laboratories at the 
Naval Postgraduate School. The various short programs developed in 
the process can be adapted easily to other geometries and to different 
applications. Extensive use has been made of the Hewlett-Packard 
Model 9830A Calculator. This type of programmable calculator allows 
performance and design calculations to be programmed very quickly by 
an individual investigator or engineer. A collection of these mini- 
programs will become the tools for a complete design analysis, and 
it is therefore advisable to document the approach as it evolves. The 
programs themselves are not listed here. It is the author's opinion 
that the programming is simple and best carried out by the user. 

However, the programs reported are identified by log number and can 
be listed or taped if requested. 

The style of the analysis results from the repeated use of the 
"total flow function" which is reintroduced in Appendix A. In the 
analysis of internal compressible flows, particularly if they are 
adiabatic, the flow function provides the means of replacing density 
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in the equations with a well-understood function of the velocity. 

With a technique for obtaining velocity from the flow function worked 
out, as described in Appendix A, the treatment of adiabatic internal flows 
as being compressible is little more elaborate than the treatment as 
incompressible, provided that a programmable calculator is available. 
Furthermore, the extreme precision of such calculators allows the use 
of the compressible flow formulation to very small velocities. 

In Section 2 of this report, results are given for a calculation 
of inviscid, irrotational but compressible flow streamlines at the 
rotor face of the TRANSX compressor. Program IPP used in these cal- 
culations was written for the IBM360 computer. When reprogrammed 
for the HP 983 OA it will be reported separately. In Section 2, the 
computed distributions of velocity and flow angle are well represented 
by a particular curve fit described in Appendix B. The curve fitting 
procedure provides the exact value at the outer wall, control of the 
slope at the wall and an exact fit through three further points. The 
technique is useful when data is known accurately at only a few 
points. In Section 3 the effect of a boundary layer on the outer 
wall is represented by a blockage factor and related to the boundary 
layer displacement thickness. In Appendix C the blockage factor for 
the TRANSX compressor is evaluated from preliminary test results. 

The analytical representation of the compressible flow in radial 
planes obtained in Sections 2 and 3 allows conditions radially along 
and relative to the rotor blading to be related to the total flow 
rate and compressor speed. In Section 4 the flow rate is related to 
the flow angle at the rotor tip for a given speed. The results 
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graphed by the HP 983 OA in Figure 4 are used to show where the compressor 
is operating in relation to design incidence conditions. In Section 
5, the distribution of the relative Mach number is calculated for a 
given flow angle at the tip and a given speed. The results graphed 
in Figure 5 are needed to guide the interpretation of flow field 
measurements as the compressor speed is increased. 

2. UWISCID, IRROTATIONAL COMPRESSIBLE FLOW STREAMLINES 

Compressible flow profiles were calculated using an IBM 360 computer 
program (IPP)^ which assumes that stream surfaces are axisymmetric and 
locally parabolic in a radial plane. The position and slope of the 
streamlines given by Program IPP agreed closely with results of a flow 
net construction for incompressible flow. [The small variation from 
hub to tip in the velocity ( 0.25 °lo) and therefore also in the density 
accounts for the agreement] . The calculation was done for the design 
flow rate corresponding to an axial Mach number at the tip of 0.646. 

It was therefore assumed that the pattern of streamlines was inde- 
pendent of the flow rate. 

Referring to Figure 1, continuity at station 1 gives 

w = f pVcos\ 2nrdr (l) 

r . 

which can be written as 1 

♦ P f r ° 

W = w/(p t V t nr 0 = J [pv/p V ]. cos\2- (r/r 0 ). d (r/r a ) (2) 

i 

where subscript t denotes "stagnation" values. 

Using the definition of the "total flow function", $, given in 
Appendix A, and writing R = r/r„ , Eq. (2) becomes 
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w 



p/ §cosXd(R 2 ) 



R i 



(2a) 



which can be written as 



w = 



= i 0 cosXoj~ (§/§ D )(cosx/cosX 0 ) d(R 2 ) 



(3) 



R. 

1 



In the TRANSX compressor X„ = 0, so that from Eq. (3)> 



* 0 = w /§ 



(4) 



where 



1 



§ 



(§/§ o ) cosxd(R ) 



(4a) 



Rr 

i 

and ^ is a function only of the inlet geometry. 

Using the data from Program IPP, the distribution of the total 
flow function could be represented analytical by 



- V 1 - 1/ 1 + "i (1 - fr> 

where the values 

*\ 

= 1.79864 

n = 0.38582 ► 

= 22.3048 



(5) 



(5a) 



were found (Appendix B) to fit the data very closely. This is 
seen in Figure 2. Similarly the flow inclination was represented 
analytically by 

n 2 

X = a 2 (l - R 2 ) + m 2 (l - R 2 ) (6) 
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and the values 



a, 



2 



O .9892665 



n. 



2 



3.8167895 



(6a) 



m. 



2 



0.4714926 



were found to fit the data very closely, and this is shown in Figure 
3. The curve fitting method is described fully in Appendix B. 

Using Eq. (5) and Eq. (6), the value of ^ can be obtained by 
integration of Eq. (4a). Since Eq. ( 5 ) cannot be solved explicitly 
for §/$ 0 as a function of R, ^ was evaluated for the TRANSX compressor 
by graphically integrating the data from Program IPP. A value of ? = 
0.708355 was obtained. 
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3. BOUNDARY LAYER BLOCKAGE 



Preliminary measurements of impact pressure ahead of the rotor have 
shown that a "boundary layer" exists on the outer wall extending appro- 
ximately 0.25" into the flow. At smaller radii the stagnation pressure 
is nearly constant. In order to correct the inviscid streamline calcula 
tion for this deficit a blockage factor can be introduced which can be 
evaluated from measurements, or assumed a priori. If 5 * is the displace 



ment thickness and b*/ T o « 1 , 
w = w + o.V a 2rr r 6 * where w 
is the measured flow rate and 
w is the flow rate for inviscid 
flow. Then, using the definition 

in Eq. (2), 

- - s* 

W = W + $ 2 - 7 - 
m ° r 




OUTER MILL 






Writing kg for the area blockage factor usually defined as 



w 



m 



w 

then 



(7) 



25 * 

W = k W + $• — 

0 r o 



or . 



V 1 



o_ _ 26^ 

r 



W 



and using Eq. (4), 



s* 



k = 1 - 

*1 S r o 



(8) 



Hence to account for the viscous blockage, kg must be assigned a 
value or calculated from data using Eq. (8), and then, using Eq. (4) 
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and Eq. (7) , 



1 Wm 




(9) 



is the value of the flow function at the tip in terms of the measured 
flow rate. (The boundary layer on the hub is very thin and is ignored 
here) . 
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h. CALCULATION OF THE FLOW RATE FOR CONDITIONS GIVEN AT THE TIP 



The above Eq. (9) relates the flow function at the tip to the 
measured flow rate. It is desirable to relate the measured flow rate 
at a given rotor speed to the flow angle or incidence at the tip. The 
"design point" is taken as the condition at which the flow is at the 
design incidence angle at the tip = P-^q ) and the "design" flow 

rate is therefore a function of the rotor speed. 

If p is given, using the notation of Fig. 1, 



If 9 = T-tVTpgfj where "ref" denotes a reference temperature (taken as 




( 10 ) 



and 



V = — 

o tanp i0 



U 

o 



( 11 ) 



T re f = 517. 8°R), then 



X 



o 




( 12 ) 



where Eq. (9) and Eq. (10) have been used 



V 



ref 




(13) 



and 



V 



t 




( 14 ) 



as described in Appendix A. 



The total flow function at the tip is given by 
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( 15 ) 



§ 



o 



2 1 / 7-1 
X (1 - x 2 ) 

O O 



and. using Eq. ( 9 ) > 



W 



m 






From the definition of in Eq. (2), 



9 

w 

m 






w 



m 



( 16 ) 



Using Eq. (l6), Eq. (15) and introducing reference pressure and 
temperature , 




f r o\ 

km~/ 



p ref^ref\l44 / ^B1 



? • 



r 2-f/y- 1 

V^oj 



(17) 



In this equation Pref = (e/fi ) > where 5 = (p t /p ref ), P ref is a 
reference pressure (taken as P re f = 29.92" Hg) and a perfect gas is 
assumed. 

Since Eq. (12) relates the outer velocity X q to the referred speed, 
N/ \fj and the flow angle P-^qj Eq. (17) gives the referred flow rate in 
terms of these two variables for a given blockage factor. The HP 983 OA 
was programmed to graph the relationship given in Eq. (17). The result, 
for the blockage factor ( 0 . 99 ) measured in preliminary tests of the 
TRANSX compressor is given in Fig. 4. (The program is logged as HP- 
Sf001-15). 

It should be observed that the design flow angle at the tip for the 
TRANSX compressor is p^ = 65 °. The design speed of the compressor is 
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30,460 rpm. In Fig. 4 choking of the inlet flow is seen to occur at 
about p = 57 " at this speed. 

5. CALCULATION OF RADIAL PROFILES FOR A GIVEN FLOW RATE 

For a given measured flow rate and stagnation conditions, W can be 
evaluated from the definition in Eq. (2). is then given by Eq. (9). 

At any value of R, $/§ q is given by Eq. (5) using the method of 
iteration given in Appendix C, so that § is known. The non-dimensional 1 
velocity X is obtained by solving the flow function as described in 
Appendix A. Then the non-dimensional local peripheral speed is 



X 



u 




(18) 



and relative speed is 



X 



w 




(19) 



The local relative Mach number is then given by 



M 



W 



W 



1 



( 20 ) 



w 




and the flow angle by 




( 21 ) 



This procedure can be used to obtain profiles of relative Mach 



numbers and flow angle for measured test data. 
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Alternatively, by specifying the flow angle at the tip can be 
obtained from Eq. (12) and Eq. (15) and the radial profiles calculated 
in the same way. 

The latter method was programmed for the HP 983 OA Calculator to 
obtain distributions of relative Mach number at "design" point of the 
TRANSX compressor as a function of the referred speed. (The program 
is logged as HP-Sf001-l6) . The results are shown in Fig. 5. It can 
be seen in Fig. 5 that at the design speed of 30,460 rpm, only a small 
region of the flow near the hub has a subsonic relative Mach number. 
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Appendix A. Velocity from the Total Flow Function 



The "total flow function" is the mass flux divided by the "limiting" 
or "stagnation" mass flux. A discussion of the flow function is given 
in Appendix A of reference 2. In addition to defining stagnation 
values of pres sure, temperature and density, a "stagnation" or "limiting" 
velocity is introduced from the energy equation. Thus if 



\ ' h + I 



A(l) 



defines the total enthalpy, 




A(2) 



is the "limiting" velocity obtained in an adiabatic expansion to 
zero temperature. 

For a perfect gas, the equations relating static and stagnation 

properties become very simple if the velocity is made dimensionless 
(3) 

by dividing by V t w/ . Writing 



then 



X 



V/V t 



T / T t 

p/p t 

u/p t 



1 - X 



2 y/(y-D 

(i-x 2 ) 



2 Vfr -1 ) 

(l-X 2 ) 






A(3) 



A(4) 
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and the total flow function, § , is then given by 



? 1/ (*y-i) 

§ = P v /p t v t = x(l - X*) A(5) 

In analysing internal compressible flows, the flow function in 
this form can often be used to eliminate density from the equations. 
In many cases, the flow rate and the stagnation properties are 
known and Equation A(5) must be solved to obtain the velocity. Care 
must then be taken to recognize the subsonic and supersonic roots as 
shown in the following sketch of Eq. A(5): 



i 




Differentiating Eq. A(5), 



d$/dX = $[l/x - 2X/(y - 1)(1 - x 2 )] A(6) 



lb 



For the subsonic root the procedure is as follows: 



(i) guess X < (y-l/y+l) ; X = $ (l+§ ) is suitable 

x , 2 1 / (y ' l) 

(ii) calculate § = X (l - x ) 

(iii) test for convergence; - $ J < e? where is the given 

value 

(iv) calculate X = X + ($^ - #)/(d$/dX) using Eq. A(6) 

(v) repeat from (ii) for convergence at (iii) 

For the supersonic root choose an initial guess at (i) such that 
X > (l - iy(y + l) and use the above procedure without change. 
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Appendix B. Curve Fitting with CT 



= ae 11 + me 

Compressor blading calculations involve the selection of a number of 
stream surfaces between which equal fractions of the total flow rate 
pass. By referring properties to conditions at the outer radius, 
distributions of quantities of order unity or less must be represented. 

A function which has been found to be useful to approximate smooth 
distributions calculated only at a limited number of points is the 
following; 

<j = as 11 + me B(l) 

The shape of a = CT (e ) given by Eq. B(l) depends on the range of 
e and particularly on the value of n. Differentiating, 



d^/d* = ane 11-1 + m B (2) 

from which it is evident that there is a discontinuity in behavior at 
n = 1 for n going from values less than to greater than unity. This 
can be seen in the following sketches: 

Case 1: YL < » Case 2: YL > I 





l6 



In Case 1, the slope is infinite at e = 0 and is given by the value 
of m as e becomes large. In Case 2, the slope is given by m at e =0 
and increases as e becomes large. 

The function can be fitted to given data in a number of ways. There 
follows the analysis for fitting three arbitrary points exactly. (In 
application it is found that particular choices of data points can 
be made to improve the overall fit.) 

If subscripts 1, 2 and 3 denote three given data points, the equations 



°1 “ 84 ]_ + me ^ 



ct 2 ~ ae 2 + me 2 



0 ^ ae ^ m e ^ 



B(3) 



must be solved to give the values of a, m and n. 

Eliminating m between pairs of Eq. B(3) and 'a' from the resulting 
pair of equations, n is given by 



-1 n-1 

1 e 2 



n-1 n-1 

s 2 “ S3 



e 



e 



■^ ( 0 

•j- ( 6 2 



e-^Og) 
e ^ ) 



B(4) 



When n is known, a and m are given by 

s ® 2^0" 2 

a = n n 

e 2 e l " e l e 2 



B(5) 
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and 



m - (a i - ae i n )/ ei = - ae^)/e 2 B(6) 

Thus it is necessary first to solve Eq. B(4). Writing 

L = [LHS of Eq. A(4)] = - 1]/[1 - (^/e,,) 1 " 1 ] 

or 

L = (A 11 - 1)/(1 - B 11 ) B(7) 

where A = B = e ^/e 2 and T| = n - 1. 

Choose 6l , e 2 and such that 0 < 6l < ^ 2 < e then A < 1, B > 1 
and T| (< 0) is to be determined. 

Differentiation of Eq. B(7) gives 

dL/dn = dL/dT| = (A^ - B 71 )/(l - B 11 ) 2 B(8) 

From an inspection of Eq. B(7) and Eq. B(8); 

(i) If T] > 0 (n > 1) , L > 0 and (dL/dn) < 0 

(ii) If T| < 0 (n < l) , L > 0 and (dL/dn) > 0 

Note that as T) -* », L _ 0. Hence the situation for obtaining a 
solution to Eq. B(4) is as shown in the sketch, where 



l8 




So ( S 1 ” S -l O p ) 

R = TRHS of Eq. B(U)l = B(9) 

€ 2_ ' s 3 a ^ 3 / 



If an initial value of 7] is chosen ( = t)^) , then an improved guess 
is given by 



T) 



\ + J 





where j = + 1 when T|^ 



+ve 
-ve ' 



B(10) 



The initial guess should be made from an inspection of the curva- 
ture required to fit the data. 

If we define the difference in slopes as follows: 



C = ct 3 ~ g 2 
e 3 " e 2 

then 



C > 0 


requires 


T| > 1 


(T) > 


0) 


O 

A 

o 


requires 


T) < 1 


('ll < 


0) 



ct i 



B (ll) 
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